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Abstract 

Background: In the classical view, cell membrane proteins undergo isotropic random motion, that is a 2D 
Brownian diffusion that should result in an homogeneous distribution of concentration. It is, however, far from the 
reality: Membrane proteins can assemble into so-called microdomains (sometimes called lipid rafts) which also 
display a specific lipid composition. We propose a simple mechanism that is able to explain the colocalization of 
protein and lipid rafts. 

Results: Using very simple mathematical models and particle simulations, we show that a variation of membrane 
viscosity directly leads to variation of the local concentration of diffusive particles. Since specific lipid phases in the 
membrane can account for diffusion variation, we show that, in such a situation, the freely diffusing proteins (or 
any other component) still undergo a Brownian motion but concentrate in areas of lower diffusion. The amount of 
this so-called overconcentration at equilibrium issimply related to the ratio of diffusion coefficients between zones 
of high and low diffusion. Expanding the model to include particle interaction, we show that inhomogeneous 
diffusion can impact particles clusterization as well. The clusters of particles were more numerous and appear for a 
lower value of interaction strength in the zones of low diffusion compared to zones of high diffusion. 

Conclusion: Provided we assume stable viscosity heterogeneity in the membrane, our model propose a simple 
mechanism to explain particle concentration heterogeneity. It has also a non-trivial impact on density of particles 
when interaction is added. This could potentially have an impact on membrane chemical reactions and 
oligomerization. 
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Background 

In the classical fluid-mosaic model of membrane, mem- 
brane components undergo isotropic random motion akin 
to Brownian motion [1,2]. In this model, since specific 
interactions between individual molecules are not consid- 
ered, resulting equilibrium distribution of components is 
homogeneous. Thus, the so-called Singer-Nicolson model 
implicitly assumes that no lateral heterogeneities within 
the membrane plane exist. 

Recently, this picture has considerably evolved towards 
a non-homogeneous distribution of cell-membranes 
components [3-5]. Experimental evidence has accumu- 
lated that cell membranes, particularly cell surface mem- 
branes, are indeed laterally heterogeneous on scales that 
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range from tens of nanometers to a few microns. These 
heterogeneities are commonly referred to as "micro- 
domains" to contrast them with the membrane macro- 
domains, the functionally differentiated surfaces of 
epithelial and other morphologically polarized cells. 
More evidence points towards the fact that some of 
these domains are enriched in various lipids such as 
cholesterol. These are sometimes called lipid rafts [6]. 

The organization of membranes into microdomains 
can biologically be interesting because microdomains 
could strongly affect membrane functions as interacting 
species are likely to be in higher concentrations in the 
domains. Indeed a wide collection ofmembrane proteins 
involved in such processes have been shown to coloca- 
lize with rafts. They are thought to play an important 
role in various cellular processes such as trafficking and 
signaling to cite but a few [7,8]. In addition, pertur- 
bations (like disruption by cholesterol removal) have 
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important effects on cell responses compared to control 
[9,10]. Surprisingly some recent works even suggest that 
the non-raftparts of the membrane are protein free [11]. 
It has been reported that membrane proteins with at 
least one transmembrane domain or with a hydrophobic 
modification are enriched in lipid rafts [12]. 

It is not a surprise that various models have emerged to 
explain the existence and stability of membrane heterogen- 
eity [13]. In particular, membrane protein diffusion alter- 
ation has been extensively studied in the context e.g 
oftrapping systems [14-16], specific interactions [17-19] or 
crowding [20,21]. Membrane-skeleton fence or membrane- 
skeleton corralling models have also been proposed. Trans- 
membrane proteins protrude into the cytoplasm and, in 
this model, their cytoplasmic domains collide with the 
membrane skeleton, inducing temporary confinement 
of the transmembrane proteins within the membrane- 
skeleton mesh [22-24]. The transmembrane proteins then 
hop to an adjacent compartment. This model was sup- 
ported by Monte Carlo simulation results [25]. 

The latter model therefore rests on the assumption that 
rafts are well defined stable structure with mixed effects 
on diffusion. Indeed constrained diffusion within bounds 
imply an increase in concentration only temporarily since 
outbound protein have low probabilities to come back. In 
addition several other works suggest a more complicated 
picture [26]. Rafts (domains) are not well defined entities. 
Several measurements seem to indicate that some proteins 
are stably associated over long but finite period of time 
with discrete domains. These domains can either diffuse 
across the cell surface [27] or are immobile, such as cell 
surface caveolae [28]. Contrary to this view is the dynamic 
assembly. In this model proteins are transiently occupying 
raft domains and undergo diffusion both inside and out- 
side the raft [29-31]. 

Being bound to a domain affects the protein diffusivity 
[32]. The same applies for lipids as their lateral mobility in 
a liquid-ordered (raft) phase is slower than in a liquid- 
disordered phase (non-raft) [29,33,34]. Indeed, several 
studies suggest that proteins and lipids undergo con- 
strained and/or slowed diffusion within rafts [30,31]. How- 
ever, individual raft proteins do not appear to undergo 
correlated diffusion with one another [35]. 

We propose in this article a simple mechanism which 
leads to the formation of protein/lipid enriched microdo- 
mains. This mechanism is based on non-homogeneous 
diffusion (NHD). Indeed, many structural constituents of 
the membrane can alter its viscous (or diffusion) proper- 
ties. As such, a membrane cannot be approximated any- 
more by a simple 2D manifold with a constant diffusion 
but rather should be described by a diffusion profile: a 
space (position) dependent diffusion function. The mem- 
brane composition is therefore by itself a source of hetero- 
geneity in the displacement of trafficking proteins. 



We will assume throughout this paper that the cell mem- 
brane has a non-spatially constant diffusion tensor that is 
temporally stable within the timeframe of the diffusion. 
These assumptions allow a greater generality since one 
never needs to assume any structurally stable component 
to define a domain. However we will not address the under- 
lying mechanism behind such diffusion variation and sim- 
ply assume it exists. Recently, several works have provided 
plausible mechanisms for viscosity alterations [36-39]. 

By solving the corresponding equation of motion, we 
show that the diffusion profile relates simply to the equi- 
librium concentration profile. In the case of punctual 
particles without interaction we show that this relation 
is extremely simple and we give a closed form of the 
equilibrium. In some cases we are able to provide a 
closed form for the whole solution and derive a FRAP 
(Fluorescence Recovery After Photobleaching) estima- 
tion of the diffusion in such conditions. 

We expand the model by performing simulations with 
non-punctual interacting particles. We show that the 
classical phase transition observed for this kind of sys- 
tem is altered by non-homogeneous diffusion - namely 
resulting in a shift in the transition diagram. This shift 
allows for extremely high concentration in the slow 
zones - much higher than for non-homogeneous diffu- 
sion or interaction separately - as well as higher cluster- 
ing properties. 

Methods 

Continuous model 

In this first model, membrane particles are independent, 
punctual and subjected to Brownian motion. The mem- 
brane is either a one or two dimensional manifold with 
periodic boundary conditions. This imposes periodic 
boundary conditions on the equations of motion. 

For simplicity's sake, we first derive the model and equa- 
tions for the ID case (equations for the 2D case are pro- 
vided afterward using another method). Therefore, we 
consider the membrane as a segment of length 2 L cen- 
tered on zero ([-L; L]); when a particle crosses one of the 
boundary it will appear at the other side). We assume that 
these particles undergo a classical Brownian motion in the 
position space (i.e. overdamped regime). To describe 
membrane heterogeneity, the diffusion coefficient depends 
on the position of the particle on the membrane. This so- 
called multiplicative noise can be efficiently modeled, in 
the ID case, with the Stratonovich formalism yielding the 
following stochastic differential equation 

dX t = V(X t )odZ (1) 

Here X t describes the position of the molecule at time 
t in the segment [-L; L], dZ/dt is the classical Brownian 
noise (with zero mean and unit variance) and V a non- 
negative 2L-periodic function. Here V is thespace- 
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dependent standard deviation of the noise. Then, the ac- 
tual diffusion (denoted as V) is equal to D{x) = V 2 (x). 
We will suppose that V is continuously differentiable 
and > 0 everywhere. 

By setting p(x,t) the probability density function of a 
particle, we obtain the associated Fokker-Planck equa- 
tion [40] derived from Eq. 1: 



dp(x.t) 1 d f, N d r _ N , 



(2) 



Looking for continuous and differentiable solutions, 
boundary conditions emerge naturally from the 2 L peri- 
odicity. Namely, p(L,t) = p(-L,t) and ^-\ L = ^r\_ L 
for all t > 0. Moreover, fluxes of particles are opposed on 
each side of the border. Writing this condition leads to 
i[T>(x)p(x,t)]\_ L = -^[V(x)p(x,t)]\ L . 

The solution of Eq. 2 can be found at equilibrium, 
yielding the constant flux /: 



J = V(x)^ x [D(x)p(x)} 



(3) 



The boundary conditions impose J(-L) = -}{L) =/= 0 
and consequently 



p(x) 



n 

V(x) 



(4) 



with Q = f^ L So, up to a normalization con- 

stant, equilibrium density is the inverse of the square 
root of the diffusion coefficient (or equal to the square 
root of viscosity). 

This result remains the same for 2D (or any higher di- 
mension): Assuming a non-homogeneous brownian mo- 
tion in higher dimension via 



dX t 
dY t 



T>(x,y) 0 
0 T>(x,y) 



dZ 2 



(5) 



in the square [-L; L] 2 and dZ/dt (i= 1,2) being the clas- 
sical Brownian noise (with zero mean and unit variance). 
The Fokker-Planck equation for the density function 
p(x,y) will be 



dp(x,y, t) 
dt 



div(VV(Vp)) 



(6) 



with the boundary conditions 

d 

— [D(x,y)p(x,y, t)}^ 
= ^ x [D(x,y)p(x,y,m ±L = 0 



(7) 



for all t > 0. To solve Eq. 6, we multiply by Vp and inte- 
grate with respect to both x and y to obtain 



mil v P 2dxd J =11 T>pdiv(VV(Vp))dxdy since §- t p 2 = 
\p% An integration by parts yields J j Vpdiv x 
(VV(Vp))dxdy=f [ _ m2 V 2 pV(Vp) - fJV(Vp)VV(Vp) 

dxdy . The first integral on the right hand side of the 
previous equation is taken over the square [-L; L] 2 and 
vanishes due to toric boundary conditions (Eq. 7), 
yielding 



ld_ 
2di 



(8) 



Since V is > 0 everywhere the left hand side of Eq. 8 is 
a decreasing function and its derivative must reach zero 
at equilibrium. 

This equilibrium therefore verifies 



J J V\\7(Vp)\ 



dxdy = 0 



meaning the gradient of Vp is zero everywhere. That is 
D 



p(x,y) 



T>(x,y) 



(9) 



The same result holds for any dimension: any diffusion 
over a non-constant diffusion profile yields, at equilib- 
rium, to a non constant concentration namely its in- 
verse. Note that as it should be, in the case of 
homogeneous diffusion (V(x) = y/Do) the equilibrium 
density (p(x) = (2Z,) -1 ) is independent of the diffusion 
coefficient D 0 . 

As a consequence zones of slow diffusion (or high vis- 
cosity) will tend to gather more particles at equilibrium 
than any faster zones. Note that we used the Stratono- 
vich formalism to describe non-homogeneous brownian 
motion of particles. This formalism make assumptions 
on the diffusion at the microscopic scale [41]. With 
other assumptions, we could have derived the same 
equations using the Ito formalism to obtain a slightly 
different result. 

Briefly: the associated Fokker-Planck would have been: 



dp(x,y, t) 
dt 



■ div(V(V 2 p)) 



yielding an equilibrium distribution 



p(x,y) 



V l {x,y) 



(10) 



(11) 



The main results hold in a stronger way (due to the 
square) in any dimensions. In the continuous model, all 
simulations described below have also been done using 
the Ito formalism for the random motion scheme. All 
the results at equilibrium described below hold qualita- 
tively in the Ito formalism but in a stronger way 
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quantitatively. Thus, to stress our point, we will focus on 
the weaker' form of our results and remain in Stratono- 
vich formalism throughout the rem' aining of the paper. 

The previous results hold for equilibrium distribution 
only. General transient solutions can be obtained for 
the one dimensional case by solving Eq. 2 entirely. 
Let x(x) = f_ L and p(x, t) = V(x)p(x, t) Eq. 2 

becomes 

dp 
dt 



l d 2 p 
2d^ 



(12) 



which is exactly the classical diffusion equation. One can 
immediately see that the same boundary conditions 
apply to p by replacing x by x. Equilibrium solution for 
classical diffusion on a circle is a constant density. So 
V(x)p(x 1 1) = p(x, t) — £1 . In addition, we obtain the 
transient solutions (on a circle) by replacing x by 
I-l W(uj m me g enera l solutions of the classical diffu- 
sion. Note that depending on the initialconditions, parti- 
cles can take a longer time to reach equilibrium in the 
slow diffusing zones. Therefore, transiently it may hap- 
pen that zones of slow diffusion will gather less particles 
than in faster zones. 

On biological membranes, transient solutions are impos- 
sible to measure. The classical solution is to estimate them 
indirectly via imaging techniques such as Fluorescence Re- 
covery After Photobleaching (FRAP). Bleaching particles 
amounts to create a new initial distribution p(x,0) =g(x). Par- 
ticles inside the bleaching are removed (actually bleached) 
while the others are untouched: g(x) = 0 for \x\ < r (for a 
bleaching beam of radius r centered on zero). In that case, 
the classical approximation of the half-time to recovery 
[42,43] is 



yr^ 
D 



(13) 



where y is a constant. Replacing r by f_ r and using 

D(r) = V 2 (r), we can account for an heterogeneous diffu- 
sion dependent half-time recovery. 



T r = y 



r du 
r V(u) 




du 



(14) 



Equation 14 indicates that the time constant of the re- 
covery is proportional to the square of the spatial aver- 
age of the inverse of the square root of diffusivity. Using 
this equation, one can, theoretically, extracts the diffu- 
sion profile from FRAP data [44], Unfortunately, for a 
general diffusion profile, no such transient solution can 
be found for higher dimensions. 



At this point, we can already make some useful 
predictions on membrane particle distribution. Indeed, 
some membrane components are known to alter 
membrane viscosity. Our model suggests that, around 
components that increase viscosity, diffusing particles 
will be in higher concentration - overconcentration - 
compared to faster zones. A good candidate could be 
cholesterol: Indeed, cholesterol is found in high 
amount on virtually all mammal cells and is known 
to rigidify membranes and cholesterol-enriched mem- 
branes display lower diffusion coefficient (by around 
one order of magnitude) [45,46]. In this case our 
model predicts that membrane cholesterol islands 
(patches) are originating membrane domains. In other 
words and without any more assumptions, stable 
cholesterol domains should imply by means of slow 
diffusion a protein-riched domain. 

In addition, our model estimates this domains to con- 
centrate more particles in proportion to the slowing. 
Indeed this overconcentration is simply proportional to 
the square root of the ratio of diffusion. For example, dif- 
fusion an order of magnitude slower should provide slow 
domains with around three times more particles. 

Conversely, when cholesterol concentration is low 
enough or high enough for the membrane to exhibit an 
approximatively constant diffusion profile, our model pre- 
dicts no domain. 

Force field interactions 

This first model is neglecting some crucial features of 
membrane diffusing particles. As such, our model 
assumes independent particles that never collide nor 
hinder each other and that cannot interact. But both 
crowding and particle interactions are known to yield in- 
homogeneous particle distribution. We expand the pre- 
vious model to take into account steric hindrance and 
possible particle interactions (e.g. protein-protein inter- 
actions) [47]. 

To study the combined effect of non-homogeneous 
diffusion and particle interaction, we used 2D particle 
simulations. As such, particle interactions must be mod- 
eled at the same scale as the diffusion itself: the meso- 
scopic scale. In order to do so, particles were subjected to 
non-homogeneous Brownian motion and short- distance 
interactions [48,49]. Assuming there are N particles, we 
introduce a new equation of motion for a particle / as: 



dyj 



(15) 



where u t j is the unitary vector of the semi-line starting 
from particle i to particle /, r t j being the distance between 
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the two particles. The interaction force / is derived from 
the 6-12 Lennard- Jones empirical potential accounting 
for van der Waals interactions: 

the force / itself being the opposite of the gradient of the 
potential /(r) = — ^r^. The two parameters G and a de- 
scribe the shape of the interaction. For a given a increas- 
ing G increases the strength of the interaction (as well as 
its range) and clusters of particles appear. The average 
size ofclusters with increasing G yields a phase transi- 
tion from N clusters of size 1 to one cluster of size 
N [50,51]. 

However, in the case of non-homogeneous diffusion, the 
phase transition will be shifted to lower interaction 
strength in the slower zones. As such, fast zones will act 
as a particle provider (well) and slow zones as particles 
sink. We expect therefore to obtain in the slow zones a 
higher over-concentration compared to the fast zones but 
also a higher degree of clustering. This over-concentration 
should be higher than both effects (interaction and non- 
homogeneous diffusion) taken separately. 

Results and discussion 

Continuous model 

We present in this section the results obtained via simu- 
lation of particles that undergo equation of motion such 
as Eq. 1. We used the Milstein numerical schema for the 
Stratonovich formalism [52], 

Xi(t k+1 ) = VAiv(Xi(t k ))z 

+ ^V(X i (t k ))V f (X i (t k ))Z 2 

Where Xi (t) is the position of particles i at time t 
(t k =(k+l)At) and V and V' are the diffusion profile 
and its gradient respectively. Various diffusion profiles 
were tested as well as different number of particles. For 
the 1 dimensional experiment 100,000 particles were 
used whereas 1,000,000 were used for the 2 dimensional 
case. We tested various time step and took the highest 
value (A£ = 0.01)) that yielded no change in equilibrium. 

Equilibrium distribution behaved as expected. The first 
set of experiments is for the one dimensional case where 
we test various diffusion profiles (continuously differenti- 
able 2 L-periodic). As an example, we present two of 
them: In one case, Figure 1, the diffusion profile is a 
Gaussian function V{x) = D 0 (l — d exp(—x 2 /co) with 
D 0 = 1, d = 0.9 and co = 0.1 (displayed in inset). The theor- 
etical distribution D/V(x) and the simulated particle dis- 
tribution are depicted by a plain line and dots respectively. 



Transients for this simulation are displayed on Figure 2 
as p(t,x) for 0< t < 10 and #G[-1:1]. The initial condition 
is p{0,x) = 6 (x - 0.5) a Dirac function centered on 0.5. 

In this case, due to the slow diffusivity the concentration 
at the peak of viscosity p(t,0.0) is the lowest value for a 
long period: transiently the slow zone is under populated. 

In addition, we simulated FRAP experiments and com- 
pared the diffusion coefficients computed through simu- 
lation to the theoretical ones as a function of the FRAP 
radius r (the FRAP experiment was centered on the 
slowest diffusion point). Figure 3 shows the results for 
three experiments, one classical homogeneous diffusion 
with V(x) = 1 and two NHD with a centered gaussian 
V(x) = 1 - de~ x2 / 0) with d = 0.9 with co = 0.1 and co = 0.01 
(shown in inset). Triangles, squares and circles are diffu- 
sion coefficient computed from the particle simulation 
for the cases V{x) — 1 , V = V(x) with 6) = 0.1 and 
co - 0.01 respectively. The lines are their theoretical 
counterparts computed using Eq. (14). 

Next, we present the results of NHD on a two- 
dimensional torus. The diffusion profile consists of 
4 randomly positioned Gaussian patches centered on 
(pc b ji) G [-L; L] 2 with l</<4 and Vi(x,y) = 
D 0 (l -dexp(-(((x-Xi) 2 + (y-yi) 2 )/co) with D 0 = l, 
d = 0.9 and co = 0.1. The overall diffusion profile is the 
average V(x, y) — \Y^=\ A(# 5 y)« Figure 4 shows the 
diffusion profile (light spectrum heightmap, Figure 4A) 
and the 2D histogram of the resulting equilibrium dis- 
tribution of particles (Figure 4B). The last panel (4C) 
compares a cross-section (# = 0.2 is constant) of theor- 
etical and simulated distributions. This cross-section 
can be visualized as a black line on the heightmap of 
Figure 4A and as a white line on Figure 4B. 

In both one and two dimensional cases, the equilibrium 
distribution of particles that undergo non-homogeneous 
diffusion matches the theoretical prediction. Slower zones 
concentrate more particles than the others. In addition 
this over-concentration isproportional to the amount of 
slowing. 

Force field interactions 

Simulations were performed for the particles with inter- 
action. To take into account steric hindrance, the nu- 
merical schema (Eq. 17) was modified by adding an 
interaction strength vector AtJ[Xi(t)). 

Collision detection was included to allow a simple 
crowding experiment (i.e G = 0 in Eq. 16). We present 
the following experiment: On a 2D square (with toric 
boundary conditions) we insert a squared slow patch (on 
the upper leftcorner) where diffusion V is r = Vq/V 
times smaller than the diffusion (D 0 ) in the remaining of 
the space: a "cholesterol" patch. Another square of the 
same surface away from the cholesterol patch is used as 
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c 
0 
Q 




Position 

Figure 1 Results for 1 d experiment. Simulation (dots) versus theoretical (plain line) distribution of particles that undergo non-homogeneous 
diffusion. The diffusion profile is a simple inverse gaussian function (in inset) with equation: D(x) = 1 - 0.9exp^-^ and a=0.1 (theoretical: 
plain line, simulations: triangles). 



control patch. Numbers of particles are counted in both 
patches at equilibrium and normalized by the total num- 
ber of particles and the size of the patch. 

Assuming cholesterol diffusion is reduced by around 
one order of magnitude [53-56], those numbers are 



displayed on Figure 5 for two conservative (i.e higher) 
values of the diffusion ratio r = V 0 /V = 2 (black square) 
and r- 1.4 (black circles). The dashed lines are the the- 
oretical predictions for the concentration in the patch in 
case of punctual particles (that is y = r). For low 




Figure 2 Transient distribution. Simulation of p{t,x) for t e [0:5] and x e [-1:1]. Setting p(0,x) = 6{x-05) as the initial condition. The diffusion 
profile is as in Figure 1 and therefore the distribution obtained near f = 5 is the same as the one at equilibrium described in Figure 1. 
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Figure 3 FRAP experiment. Results for FRAP experiments compared to theoretical prediction. Three experiments with different diffusion profiles 
are depicted. D = 1 (theoreticahplain line, simulations: triangles) and two gaussian profile D(x) = 1 - 0.9exp^-^ with a = 0.01 (middle, 
theoretical: dotted line, simulations: squares) and a = 0.1 (bottom, theoretical: dashes line, simulation: circles). 



interaction strength, the punctual particles prediction 
remains a good approximation. However when the inter- 
action increases, concentration of particles in the choles- 
terol patch increases dramatically while there is no such 
increase in the control patch. On Figure 6 screenshots of 
particles positions are displayed (r = 0.5) for G=0 (A), 
G = 500(B), G = 4000(C) and G = 10000(D) - see also sup- 
porting videos. 

As expected, the overconcentration effect due to non- 
homogeneous diffusion is dramatically increased via par- 
ticles interaction. As Figure 5 shows it, a great portion 
of the particles is concentrated in the cholesterol patch. 
In this extreme case, the patch is virtually composed of 
one big cluster of particles in interaction whereas there 
are no clusters in the remaining of the space. 

To numerically assess this situation, we compute the 
clustering coefficient as a function of the interaction 
strength: c{f)=N c /N i where N+N c is the number of 



clusters and N t is the number of particles involved. Clus- 
ters were computed by creating a graph of particles whose 
nodes are the particles themselves and link between two 
nodes exists if the distance between their centers is below 
ar 0 (a = 2.7 was chosen a bit higher than a =2.5). A value 
close to 1 indicates no clustering while a more pro- 
nounced clustering yields smaller values. Results are dis- 
played in Figure 7. The values are computed in three 
situations: in the first two, diffusion is homogeneous with 
diffusion V — 1 (squares) and V = 0.5 (circles). The 
remaining case (triangles) is the one with a cholesterol 
patch with ratio r = 2 ( V 0 = 1 as in the previous exp- 
eriments). In all cases, the clustering coefficient decreases 
with the interaction strength and thereforeclusters appear 
(one needs higher interaction strength for the case V = 1). 
Moreover, the lower the viscosity the deeper the clustering. 
However, for high enough interaction strength, the choles- 
terol patch does not behave as an homogeneous slow zone: 





Figure 4 Results for 2 d experiment. A) The diffusion profile in light spectrum (toward red is slower) for 4 randomly "cholesterol" islands. B) 
2 d histogram of simulated particles positions at equilibrium. C) Comparison with theoretical prediction for a cross section (x=0.25). This section 
is represented by a black line on the heightmap and a white line on the 2 d histogram. 

V ) 
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Interaction Strength 

Figure 5 Particles interactions results. Concentration of particles in the "cholesterol" patch for two diffusion ratios r = D 0 /D with r = 2 (black 
square) and r= 1.4 (black circle). Plain line is the control concentration with D = 1. Dashed line is the theoretical NHD concentration. 




Figure 6 Particles maps. Position map of particles for A) e = 0 and r = 2, B) e = 500 and r = 2, C) e = 4000 and D) e = 1 0000. Note that the actual 
values of £ are scaled down by 10" 12 . See also supporting Additional file 1: Video SI, Additional file 2: Video S2, Additional file 3: Video S3 and 
Additional file 4: Video S4 for re{2,10} and GG{0,1e5}. 
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Figure 7 Effects on clustering. Clustering coefficient (see text for definition) as a function of the interaction strength e for control D = 1 
(squares), slow experiment D = 0.5 (circles) and cholesterol square patch withration r = 2 (triangles). 



whereas cluster sizes seem to reach a plateau in the slow 
diffusion experiment, clusters were bigger in the choles- 
terol patch. 

To mimic cholesterol effect, we deliberately chose 
diffusion ratios r 2 = V 2 Q /V 2 = Dq/D to be conservative 
i.e lower than one order magnitude of the diffusion ratio 
(r 2 ranges from 2 to 4 instead of 10) in almost all experi- 
ments [53-56]. However, for completeness' sake, we finally 
test several values of diffusion ratios. Results are summed 



up in Figure 8 for three values of the interaction strength 
(e) and for ratios going from r = 1 (no slowing) to r= 10. 
Typical extremal dynamics are available on the supporting 
Additional file 1: Video SI, Additional file 2: Video S2, 
Additional file 3: Video S3 and Additional file 4: Video S4. 

These results first confirm the predictions for punctual 
particles. With the parameters tested and as shown on 
Figure 8, the crowding only experiment did not differ much 
from the predictions. In addition, we show that particle 




2 4 6 8 10 

Diffusion Ratio 



Figure 8 Diffusion value dependence. Concentration ratio as a function of the ratio of diffusion for three interactions strength e = 0 (triangles), 
f = 1000 (squares), and, e = 10000 (circles). Bold dashed line is the theoretical prediction of concentration for the case of no interaction (punctual 
particles) y = x. To guide the eye the line y= 1 is displayed. 
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interactionsin the form of Lennard-Jones potential can 
strongly modify both the magnitude and the nature of the 
resulting overconcentration. Around conservative values 
for diffusion, the so-called cholesterol patch can gather up 
to 9 times more particles than expected for an homoge- 
neous medium. In addition, these particles tend to be more 
in clusters than any other situation. Our results can only be 
obtained and explained by combining both NHD and par- 
ticle interactions. The cholesterol patch - the domain - 
spontaneously gathers particles and increases the stability 
of their interactions. 

Conclusion 

Cell membrane can display a wide variety of heterogeneity 
in its physical properties, viscosity and constituents com- 
position. Our model describes the possible emergence of 
what we call a minimal membrane domain based on only 
one feature of the membrane, namely its local viscosity. We 
show that, everything else being equal, domain with high 
viscosity tends to gather statistically a larger proportion of 
the diffusive particles. Particles are not trapped within these 
domains but merely tend to spend more time in them. In 
addition, we show that the ratio of concentration of parti- 
cles in and out domain is proportional to the square root of 
the inverse ratio of their local diffusion. As such, protein 
free membrane zones are seen as a product of high diffusiv- 
ity and not from forbidding constraints - contrary to the 
hypothesis developed in [11] - and in accordance with the 
results found in [26]. Provided we assume that viscosity is 
indeed non constant for relevant time scales, this provides 
an extremely simple and parsimonious explanation for 
protein-enriched membrane domains. 

Moreover, while this model is fairly general, we 
emphasize the particular role of cholesterol. Cholesterol 
is ubiquitous in cell membranes and its influence 
on membrane diffusion has been well documented. 
Therefore, it is not surprising, from our point of view, 
that cholesterol has been found in higher quantity in vir- 
tually all microdomains exhibited so far. For many 
authors, its presence even characteries membrane 
domains [3,11]. Our model shows that gradient in chol- 
esterol concentration leads naturally to protein domains. 
This property depends only on the relative diffusion 
between cholesterol-rich zones and the remaining of the 
membrane. 

Looking for experimental results, authors have 
reported the impact of temperature and cholesterol de- 
pletion on membrane heterogeneities. Cholesterol re- 
moval should decrease the colocalization of membrane 
components as in [57,58]. In addition, cooling decrease 
the diffusion ratio between membrane with and without 
cholesterol. For low temperature it even flips: our model 
predicts no colocalization or domains in the case of low 
temperature as it has been reportedin [58,59]. 



Important diffusion variation leads not only to a quan- 
titative particles concentration variation but also to a 
qualitative variation. Protein interactions alter signifi- 
cantly the particles concentration landscape when com- 
bined with NHD. We show that cluster of particles are 
more stable in slower zones and appear for lower values 
of interactions. Our model predicts that domains of slow 
diffusion will alter affinity interactions. That is, for ex- 
ample oligomerization will happen preferentially inside 
these domains. This latter prediction has been observed 
on model membrane where modifications of the choles- 
terol content trigger oligomerization [60]. Thus, chem- 
ical reactions will have their equilibrium shifted in these 
zones, revealing a potential function for membrane 
microdomains. 

Additional files 



Additional file 1: Particle evolution for r=2 and no interactions 

(e = 0). Evolution under non-homogeneous diffusion of 1,000 particles 
(depicted in green) for 1000 time step. A square slow patch on the upper 
left corner has a diffusion ratio of 2 and particles interaction is limited to 
collision. This is the typical NHD simulations when, at equilibrium, there 
will be twice more particles in the slow patch. 

Additional file 2: Particle evolution for r=2 and no interactions 
(€ = 10000). Evolution under non-homogeneous diffusion of 1,000 
particles (depicted in green) for 1000 time step. A square slow patch on 
the upper left corner has a diffusion ratio of 2 and particles interaction is 
set to an intermediate level e = 10000. For this important level of 
interactions, clusters appear everywhere but they are always more 
pronounced in the slow patch. 

Additional file 3: Particle evolution for r= 10 and no interactions 

(e = 0). Evolution under non-homogeneous diffusion of 1,000 particles 
(depicted in green) for 500 time step. A square slow patch on the upper left 
corner has a diffusion ratio of 10 and particles interaction is limited to 
collision. The important amount of slowing allows for an important over- 
concentration (that will ultimately be 10 fold). 

Additional file 4: Particle evolution for r= 10 and no interactions 
(€ = 10000). Evolution under non-homogeneous diffusion of 1,000 
particles (depicted in green) for 500 time step. A square slow patch on 
the upper left corner has a diffusion ratio of 10 and particles interaction 
strength is set to 10000. In this last experiment,the over-concentration 
appears first on the border of the patch. The difference in the phase 
transition is extremely more pronounced in the slow patch. 
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